Preview onto version 0.3.4

Here you can do some experiments with Rascal. The computation time is restricted to about 10 secs to prevent high load on server. Have a look at the documentation or at the examples. If you see wrong plots reload may help. You job is cancelled as soon someone else accesses this page, so consider a local installation.

Rascal-toolbox-0.3.4 (C) 2001,2002 Sebastian Ritterbusch

Enter 'help' for help or 'licence' for information about licence and warranty.

 "rascal.rc loaded, predefined e, i, pi"
>"Linear interval-vector iteration";
>
>Q2(A)=[A(1,1) A(2,1);A(2,1) -A(1,1)]/sqrt(sqr(A(1,1))+sqr(A(2,1)));
>
>S2=1/sqrt([2,2]);
>Precision(1,-4);
>A=[S2 S2;-2 S2]
 [[0.7071,0.7072] [0.7071,0.7072];-2 [0.7071,0.7072]]
>
>X0=Yb=Xb=Xa=[[1,2];-1];
>B=1;
>i=0;
>
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [1 [0.7071;2] [0.7071;2] [0.7071;2]]
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [2 [1.9142;2.8284] [0.9142;2.8284] [0.9142;2.8284]]
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [3 [3.3536;5.8284] [2.6464;0.1716] [2.6464;0.1716]]
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [4 [6.4926;10.8284] [1.9926;5.1716] [1.9926;5.1716]]
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [5 [12.2478;2.0642E+001] [2.2478;7.6421] [2.2478;7.6421]]
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [6 [2.3257E+001;3.9092E+001] [6.9933;0.9081] [6.9933;0.9081]]
>Bn=Q2(A*B);Yb=inv([1,1]*Bn)*A*B*Yb;B=Bn;[i++ diam(Xa=A*Xa) diam(Xb=B*Yb) diam((A^i)*X0)]
 [7 [4.4087E+001;7.4156E+001] [5.5871;13.3444] [5.5871;13.3444]]
>

Questions, Problems, Bugs? Please use this form.